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^ ' Abstract 

The stochastic hmit for the system of spins interacting with a boson field is investigated. 
^ ■ In the finite volume an application of the stochastic golden rule shows that in the limit the 

r — , dynamics of a quantum system is described by a quantum white noise equation that after 

I taking of normal order is equivalent to quantum stochastic differential equation (QSDE). For 

the quantum Langevin equation the dynamics is well defined and is a quantum fiow on the 
infinite lattice system. 



1 Introduction 



P5 ■ Starting with the work of Glauber |]TU| the dynamics of infinite classical lattice systems has been 



Q^! class of classical Markov semigroups (cf. for a general survey and for further references). 



considered by many authors and has led to study the ergodic and equlibrium properties of a new 



Quantum analogues of these semigroups have also been considered by several authors (e.g. 

11121 , ||15| , ||16|| , . . . ). However the problem of deriving these Markovian semigroups, and more 



■ generally the stochastic flows, as limits of Hamiltonian systems, was open both in the classical 



and in the quantum case. 

On the other hand the stochastic golden rules, which arise in the stochastic limit of quantum 
theory as natural generalizations of Fermi golden rule 0, provide a natural tool to associate 
a stochastic flow, driven by a white noise equation, to any discrete system interacting with a 
quantum field. Moreover another quite general result of the stochastic limit is that the Markov 
semigroup, canonically associated to the flow, leaves invariant the abelian subalgebra generated by 
the spectral projections of the Hamiltonian of the discrete system so that the associated Markov 
process describes the jumps among these energy levels. 

Since a quantum spin system in finite volume is obviously a discrete system, the above results 
naturally suggest the conjecture that, by coupling such a system to a quantum field via a suitable 
interaction, applying the stochastic golden rule and taking the thermodynamic limit, one might 
obtain a class of quantum flows which, when restricted to an appropriate abelian subalgebra, gives 
rise to the classical interacting particle systems studied in classical statistical mechanics. 

In the present paper we will prove this conjecture for general finite range interactions thus 
extending a previous result obtained by the authors in the case of the Ising model . The model 
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we consider is essentially the same as the one considered by Martin and Buffet apart from 
the minor difference that these authors consider fermion reservoir and we a boson one. Since the 
stochastic golden rule holds also in the fermion case there is no difficulty in extending the 

present results to the fermion case. We investigate a quantum system with a Hamiltonian of the 
form 

H = Ho + \Hj 

where Hq is called the free Hamiltonian and Hi the interaction Hamiltonian. The general idea 
of the stochastic limit (see P]) is to make the time rescaling t t/\^ in the solution of the 
Schrodinger (or Heisenberg) equation in interaction picture t/j'*''' = ^^^Ho^-itu ^ associated to the 
Hamiltonian H, i.e. 



with Hj{t) = e**^''if/e~**''^°. This gives the rescaled equation 



and one wants to study the limits, in a topology to be specified, 

ltot/™=C/,; limiff,(±)=ff. (2) 

After the rescaling t t/\^ we consider the simultaneous limit A — > 0, t — > oo under the 
condition that tends to a constant (interpreted as a new slow time scale). This limit captures 
the dominating contributions to the dynamics, in a regime of long times and small coupling, arising 
from the cumulative effects, on a large time scale, of small interactions (A — 0). The physical idea 
is that, looked from the slow time scale of the atom, the field looks like a very chaotic object: a 
quantum white noise, i.e. a (5-correlated (in time) quantum field 6*(t, k), b{t, k) also called a master 
field. 



2 Physical model 

In the present paper we consider the open Ising model, describing the interaction of a system S of 
spins (or, more generally, two-level systems) with a reservoir, represented by a bosonic quantum 
field. The total Hamiltonian is 

H = Ho + XHi = Hs + Hr + XHi 

where is the free Hamiltonian of a bosonic reservoir R: 

Hr = J uj{k)a*{k)a{k) dk 

acting in the representation space JF corresponding to the bosonic equilibrium state at temperature 
and chemical potential n and the assumptions we need on uj{k) are given at the beginning of 
Section HI below. Thus the reservoir state is Gaussian with mean zero and correlations given by 

(a-(fc)a(fcO) = ^,(,(,)'^,_^ ^(A;-fcO 
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The spin variables are labeled by the lattice Z*^, and, for each finite subset A C Z'^, the system 
Hilbert space is 

and the system Hamiltonian has the form 

where cr^, cr^, are Pauh matrices (r e A) at the r-th site in the tensor product. For any r, 
s e A 

Jrs — J sr ^ ^1 J'rr 

In the present paper we consider a system Hamiltonian that describes the interaction of spin 
with a finite number of other spins (finite range potential). The simplest example is the nearest 
neighbor interaction (Ising model), discussed in Section V below. The interaction Hamiltonian 
Hi (acting mi-is® T) has the form 

Hi^Y. < ® H9r), Hg) = Mg) + A*(g), Mg) = / dkg{k)a{k), 

where ip is B, field operator, A{g) is a smeared quantum field with cutoff function (form factor) 
g{k). To perform the construction of the stochastic limit one needs to calculate the free evolution 
of the interaction Hamiltonian: Hi{t) = e'^^^°Hie~'^^^° . We use the formula 

£,. = ±1 

where \sr) is the eigenvector of corresponding to the egenvalue Er to rewrite the interaction 
Hamiltonian in the form 

Hi = i: E {er,CT:iir)\er){l^r\®(A(g) + A*(g))^(DA + Dl)(A{g) + A*{g)) 

where 
Since 

[HsAlr){-lr\] = A{r)\lr){-lr\ , A(r) = - ^ J,,a,^ 

s 

the free evolution for Dj. has the form 

The sum over s here is finite because the spin at r interacts only with a finite number of spins 
(only a finite number of Jfs for fixed r is nonzero). 

Using the formula for the free evolution of bosonic fields 

e'*^«a(A;)e-'*^« = e-''''^^^ a{k) 
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we get for the free evolution of the interaction Hamiltonian: 

Hi{t) = E / dkgik) (D,e-'*("('=)~^('^))a(A;) + D;e-''^''^''^+^^'^^a{k)) + h.c. (3) 

and, for the free evolution of Dr'. 

= n E e-'''-''\6s){6s\ (4) 

seS(r) £s=±l 

where S(r) denotes the (finite) set of spins interacting with the spin r (so r 7^ 5](r)). When 
no confusion is possible we use the same symbol S(r) to denote also the set of configurations of 
correponding spins: cr(r) = {es}(s G S(r)). 
We denote 

E{r,a{r))= E JrsSs = Hs{a', = -l,a{r)) - Hsia^: = -l,a{r)) (5) 

the energy difference between two configurations with all spins fixed, and given by a"(r), with the 
exception of Er that changes value from —1 to 1. 
With these notations formula (|^) takes the form 

Hj{t) = J2 fdkg{k)(i: Dr n ks)(^s|e-^*("('')+^('^'^('-)))a(A;) + 

reA VGS(r) seS{r) 

+ E n ks)(£s|e-^*("(')-''(^'"('^)))a(A;)] +/I.C. = 

aGS(r) seE{r) ^ 

= E E /rfA;^(fc)(a,.we-^*(-('=)+^(^''^(^)))a(fc) + 

rGAo-(r)GS(r) ^ 

+G;,(,)e~'*("('=)-^('^''^('')»a(A;)) + /i.e. = 
= E E /f^A;^(A;)F,%(,)e-'*("('=)-^('-''^('^«)a(fc) + /i.c. (6) 

rGA(T(r)6S(r) 

where we denote 

Gr,a{r) = Dr ]^ \es){es\, Fr^ai^r) = G*r-a(r) ^ Gr,a{T) (7) 

To prove (P) we used a suitable rearranging of the indices cr(r) and the property 

E{r, —a{r)) = —E{r, 0"(r)) 
where — cr(r) is the configuration of spins with all spins opposite to the spins in cr(r), i.e. es(— cr(r)) = 

The operator Fr^a{r) flips the spin at site r (or kills a spin configuration). 
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3 The stochastic limit of the model 

In this section we will denote i? G 5 the pair of indices {r,a{r)) where r takes values in A. In 
these notations the free evolution of the interaction Hamiltonian takes the form 

Hi{t) = E / c^A;^(A;)F*e-'*("('=)-^(^))a(A;) + h.c. (8) 

In the stochastic limit the field Hj{t) gives rise to a family of quantum white noises, or master 
fields. To investigate these noises, let us suppose the following: 

1) iu{k) > 0, VA;; 

2) The d — 1-dimensional Lebesgue measure of the surface {k : uj{k) = 0} is equal to zero (so 
that 5{uj{k)) = 0) (for example uj{k) = k'^ + m with m > 0). 

Now let us investigate the limit of Hj{t/X'^) using one of the the basic formulae of the stochastic 
limit: 

lim 1 exp (^/(^)) = 2n6mfik)) (9) 

Since the term 6{f{k)) in (^ is not identically equal to zero only if f{k) = for some k in a set of 
nonzero d— 1-dimensional Lebesgue measure (in our case f{k) = uj{k) — E{R)), in the stochastic 
limit of the rescaled interaction Hamiltonian (^ only the terms with i? G 5+ will survive. Here iS+ 
(respectively iS_) denotes the spin configurations (e^) whose mean energy difference I]ses{r) Jrs^s is 
strictly positive (respectively negative). We will call such spin configurations positive (respectively 
negative). 

The rescaled interaction in (|^) is expressed in terms of the rescaled creation and annihilation 
operators ax^ji{t,k) = je~^^^'^^''^~^^^^^a{k). After the stochastic limit every rescaled annihilation 
operator correspondent to positive spin configuration generates one non-trivial quantum white 
noise 

bn{t,k) = Uma,,n{t,k) = lim ^ e-^^(-('=)-^(^))a(A;) 

with the relations 

[bnit, fc), b*j,{t', k')] = lim[aA,ij(t, fc), al,^(t', k')] = 
= lim 4 e-''#("(')~^(^»5(A; - k') = 2^6(1 - t')6{uj{k) - E{R))6{k - k') (10) 

[hR{t,k),hl,{t',k')]=Q 

by d^). Moreover, configurations, corresponding to different values of the energy difference E{R), 
are independent. For generic interactions ii R' then corresponding energy differences E[R) ^ 
E[R') and correponding white noises are independent. 

The stochastic limit of the interaction Hamiltonian is therefore equal to 

h{t) = lim J Hj = J dkg{k) Y: nMt, k) + h.c. (11) 

Physically such a form of the evolved interaction Hamiltonian can be explained as follows: the 
operator Fr, where R E S^, decreases the energy of the spin configuration. Therefore the vertex 
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g/3(a;(fc)-^j) _ I 
1 



F^a{k) describes the absorption, by the spin configuration R, of an energy quantum of the field 
of momentum k. 

The state on the master field (white noise) bait, k), corresponding to the equihbrium state of 
the field, is the mean zero gaussian state with correlations: 

{b*j^it, k)bR{t', k')) = 2n6{t - t')5{uj{k) - E{R))5{k - k') 

{bn{t, k)bUt', k')) = 2n6it - mu{k) - E{R))6{k - k').^ _ 

and vanishes for noises corresponding to different positive spin configurations. 

Now let us investigate the evolution equation in interaction picture for our model. According to 
the general scheme of the stochastic limit, up to possible divergences (due to the thermodynamic 
limit) that we will discuss later, we get the (singular) white noise equation 

j^Ut = -ih{t)Ut (12) 

whose normally ordered form is the quantum stochastic differential equation ^ 

dUt = {-idH{t) - Gdt) Ut (13) 
where h{t) is the white noise (0) given by the stochastic limit of the interaction Hamiltonian and 

dHit) = Y: {F*RdBnit) + F^rfi?,V(,)(t)) (14) 

RGS+ 

f pt+dt 

dBnit) = J dkg{k) J bR{T,k)dT (15) 

According to the stochastic golden rule, ( |T^ ) is obtained as follows: the first term in (|T^ is 
just the limit of the iterated series solution for (|l]) 



T 



1 rt+dt 



The second term Gdt, called the drift, is equal to the limit of the expectation value in the reservoir 
state of the second term in the iterated series solution for (|l|) 



I i-t+dt 

hm — / 

A^O A2 Jt 



Making in this formula the change of variables r = ^2 — we get 

1 rt+dt 1-0 / f 1 \ / f 1 r \ 

Computing the (Gaussian) expectation value and using formula (H) and the fact that the limits 
of oscillating factors of the form limA^o vanish unless the constant c is equal to zero, we see 
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that we can have non-zero hmit only when all oscillating factors of a kind (with ti) in 
cancel. In conclusion we get 



Res 



g/3(<^(fc)-/.) 

and therefore, from the formula 



e'*" = = 7i5{u) - iP.P.- (17) 

-oo LJ — ijj 



we get the following expression for the drift G: 



J2 I dk \g{k)\'( ^^^^ _ ^^^^ ^_^(^(,)_^) + 



Res 



^FrF*^ 1 



uo{k) - E{R) + iO e/^C^W-A*) - 1, 

= E {FhFR{9\9)R + FrF*M9)1) (18) 
Res 

We see from formula (|18D that our model also exhibits the Cheshire cat effect discovered in P] in 
a simpler model, that is: the part of the terms in the drift (^) with i? G iS+ comes from quantum 
white noises (or, correspondingly, from stochastic differentials) and describes the self-interaction 
of such noises, but products of pairs of operators F|, with the index R corresponding to negative 
spin configuration describes self-interaction of virtual noises (corresponding to vertices without 
conservation of energy) . 

The drift term contains sums over R which are divergent for large A. Therefore in the 
Schrodinger picture we will get a divergence in the thermodynamic limit. In the following section 
we will consider the evolution in the Heisenberg picture (the Langevin equation) and we will show 
that, in this context, no divergence arises. 



4 The Langevin equation 

Now we will find the Langevin equation, which is the limit of the Heisenberg evolution, in in- 
teraction representation, of a system observable. Let X be a local operator acting on the spin 
degrees of freedom, i.e. one acting only on a finite number of spins. The Langevin equation is the 
equation satisfied by the stochastic flow jt, defined by: 

3t (X) = U:XUt 

where Ut satisfies equation (|^) in the previous section, i.e. 

dUt = {-idH{t) - Gdt) Ut (19) 
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G = E {nFRi9\9)R + FRF^i9\9)R) (20) 
Res 

To derive the Langevin equation we consider 

djt {X) = jt+dt {X) - jt (X) = du;xut + u^xdUt + du;xdUt (21) 

The only nonvanishing products for quantum stochastic differentials are 

dBR{t)dB},{t) = 2Re {g\g)Rdt, dBl{t)dBR{t) = 2Re {gWndt. (22) 
Combining terms in (^TD and using ([T9|) , (p^ , (120|) and (^) we get the Langevin equation 

djt{X) = Y.Jto e^{X)dM^{t) = E ft„,j(X)dM"^(t) + Jt o ^^o(X)dt (23) 

a n=-l,l;_R65+ 

where 

rfM-i'^(t) = rfS^(t), ^^-i,k(X) = -z[X, F^], G S+ (24) 
c?M^'^(t) = rf5^(t), 0i,^(X) = -2[X,Fr], i?G5+ (25) 

and 

^o(X) = E {<r"^ + C'^ + + O (^) = (26) 

Res 

= E (-^Im (^?|^7)i?[X, F*^Fr] + tlm ig\g)j,[X, + 
Res ^ 

+Re (^71^/)^, (2F^XF^ - {X, F*F^}) + Re ig\g)+ (2F^XF* - {X, F^F*}) 

is a quantum Markovian generator. Notice that the factors Ile{g\g)^ are > only for i? G iS+ 
and vanish for i? G 5_. 

We will also use the following notions 

<i = Re ig\g)RK,,R, «:i,«(X) = 2F^XFn - {X, (27) 

= Re ig\g)1,n_,,R, k_,,r = 2FrXF*^ - {X, FrF*^}. (28) 

Let us prove that generators above, i.e. the ^-maps satisfy some special properties which will 
be crucial to prove the existence of the solution of the flow equation. Omitting, for simplicity, the 
indices R, we will consider all the above defined operators 6o, 6i, (= 61), 6q'^^\ 9q'^^\ ki, k^i, 
acting on the Uniformly Hyperfinite (UHF, cf. 0) C*-algebra Bs generated by the identity and 
the local operators (those acting only on a finite number of spins). On this algebra the limits, as 
A —>■ Z'^, of all the structure maps are well defined and their domains contain, as a common core, 
the unital dense *-subalgebra Bq of the local operators which is also invariant under the square 
root. All our operators (as well as the operators F and F*) map Bq into itself. Moreover 6i, 
are nonsymmetric derivations and 6^'"^^^ are symmetric derivations. From the above described 
properties and from (3.2.22) of ^ it follows that ki, k^i are (symmetric) dissipations. But in 
fact Ki, have a much stronger properties, than dissipativity. We will use this property in this 
section without further comments. The crucial point is the following lemma. 
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Lemma 1. For arbitrary local X, Y, i = ±1 we have 

Ki{XY) = K,{X)Y + Xk,{Y) + 2e^,{X)e,{Y) (29) 



Proof. Let us prove (|29|) for i = — 1 (for i = 1 the proof is analogous). Formula (|28|) gives 

/t_i(T) = 2FTF* - TFF* - FF*T 

so that 

fi;_i(T) = [F, T]F* + F[r, F*] 

For T = XF we get 

K„i(Xr) = [F,XY]F* + F[XY,F*] = 
= [F, X]YF* + X[F, Y]F* + F[X, F*]Y + FX[Y, F*] = 
= [F, X]F*Y + X[F, Y]F + F[X, F*]Y + XF[Y, F*] + 2[F, X][Y, F*] = 
= {X)Y + (r ) + +2{-t) [X, F]{-t) [Y, F*] 

From (p^, ( pSf ) it follows that this is exactly (|29| ) and that proves the lemma. 
Theorem 1 . For any pair of local operators X,Y, the structure maps in the Langevin equation 



i\23i ) satisfy the equation 



9^{XY) = 9^{X)Y + X9o.{Y) + Y.c^^^ep{X)0^iY) (30) 

where the structure constants c^^ is given by the Ito table 

dM^{t)dM\t) = ^ cl^dM^it) (31) 

a 

The conjugation rules of dM"'{t) and 9a are connected in such a way that formula defines a 

*-fl0W (* O = jf-O * ). 

Proof. Follows from Lemma 1 and formulae (123), (E^), (PB|) by direct calculation. 



The existence of the infinite volume dynamics and and its approximation by finite volume 
ones, is discussed in [0, p|. 

5 The simplest case: one dimensional nearest neighbor 
interaction 

In this section in order to compare our results with the known results on Glauber dynamics 
we consider the simplest case of one dimensional translationally invariant Hamiltonian with the 
nearest neighbor interaction 

Jrs Jr+l,s+li Jrs Jr,r+1 J ^ 

In this case for every r we have 4 configurations cr(r) G S(r) of nearest neighbors of the spin at r. 
We will denote these configurations ++, H — , — h and (the first symbol is the orientation of 
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the spin on the left of r and the second — on the right). Only the configuration ++ will give a 
contribution in the stochastic limit of interaction Hamiltonian (will lie in S+(r)) and the energies 
E = E{r, a) = 2 J for different r will be equal. Then takes the form 

e,^l=++,+^-+, 

reA 

and the other operators Fj:^'^^ are defined correspondingly. In this case there is only one quantum 
white noise, denoted b{t, k), which is the stochastic limit of 

A 

and satisfies the commutation relations 

[b{t, k), b*{t', k')] = 27r(5(t - t')6{uj{k) - 2J)6{k - k') 
The Langevin equation ( ]23|) takes the form 

djt{X)= J2 3to0a{X)dM''{t) 

a=-l,0,l 



where 



dM-\t) = dB{t), 0_i(X) = ~i[X, 
dM\t) = dB*{t), 01 (X) = -t[X, 
0o(X) = (e-n en + 0«) (X) = 

E (-^Im i9\9)ls,)[X, Frrr] + .Im {g\g)l^)[X, FTFr]) + 



+Re{g\g)-^^^^ {2Fi^^>XFt^^ - {X, F^^^ F^^^^} 
+Re{g\g)l^, (2il++)XFf - {X, Fi^^^ Fi^^>} 
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